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Abstract 

We discuss the decay of B c mesons into two light mesons (tt, K^*\ p, u, (j)). All these 
decay channels come from a single type of diagram, namely tree annihilation. This al- 
lows us to derive extremely simple SU(3) relations among these processes. The size of 
annihilation contributions is an important issue in B physics, and we provide two dif- 
ferent estimates in the case of non-leptonic charmless B c decays, either a comparison 
with annihilation decays of heavy-light mesons or a perturbative model inspired by QCD 
factorisation. We finally discuss a possible search for these channels at LHCb. 



1 Introduction 

The investigation of the properties of the B c meson started in 1998 when the CDF collab- 
oration observed 20.4 events containing a B c in the channel B c — > J/iplv pQ. Since then, 
its mass and width have been measured [2], and bounds on some non-leptonic channels 
have been set {J/ip with one or three pions, D* + D°. . .). From the theoretical point of 
view, the B c meson shares many features with the better known quarkonia, with the sig- 
nificant difference that its decays are not mediated through strong interaction but weak 
interaction due to its flavour quantum numbers B — — C — ±1. Theoretical investiga- 
tions have been carried out on the properties of the B c meson, such as its lifetime, its 
decay constant, some of its form factors [3], based on OPE [U |5], potential models [6j [7], 
NRQCD and perturbative methods El EDI HQ , sum rules [121 H31 HH US] , or lattice gauge 
simulations [TUl El HE1 EE] • The properties of the B c meson will be further scrutinized by 
the LHC experiments; the high luminosity of the LHC machine opens the possibility to 
observe many B c decay channels beyond the discovery one, in particular at LHCb. 
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This article is focused on the two-body non-leptonic charmless B c decays. Indeed, the 
charmless B c decays with two light mesons (ir, , ?7 (/) ,p, cu,<f)) in the final states come 
from a single diagram: the initial b and c quarks annihilate into a charged weak boson 
that decays into a pair formed of a u and a d/s quark, which hadronise into the two light 
mesons. This picture is rather different from processes such as B c — > J/ tpn for which the 
initial c quark behaves as a spectator. The recent high-precision measurements of the 
B u ,d,s an d D s decays indicate that such annihilation processes can be significant, contrary 
to the theoretical expectation of its suppression in the heavy-quark limit. Indeed, fits of 
the data not taking into account annihilation processes are generally of poor quality. 

But the understanding of these contributions remains limited. The theoretical com- 
putation of annihilation diagrams is very difficult, so that the annihilation contributions 
are often considered as a free parameter in these decays. On the other hand, in many 
B U) d, s decays, these annihilation contributions come from several different operators (tree 
and penguin), and they interfere with many different other (non-annihilation) diagrams, 
making it difficult to obtain an accurate value of annihilation by fitting experimental data. 
For this reason, the processes such as B d — > K + K~, B d — > D~K + ,B U — > D~K°, where 
only the annihilation diagram contributes, have been intensively worked out while the 
current experimental measurements are still of limited accuracy. 

The non-leptonic charmless B c decays which we discuss in this article can have an 
important impact on this issue. We have 32 decay channels which come from annihilation 
only, as mentioned above. Moreover, these decays involve a single tree operator, which 
allows us to derive extremely simple relations among the different decay channels. Finally, 
when LHCb starts running and observes the pattern of (or at least, provides bounds for) 
the branching fractions of these decays channels, it will certainly help us to further improve 
our understanding of the annihilation contributions. 

The remainder of the article is organized as follows. In section [21 we introduce the 
decay processes which we consider in this article. In section [3j we exploit SU(3) flavour 
symmetry to derive relations among amplitudes for non-leptonic charmless B c decays. In 
section HI we discuss theoretical issues related to this annihilation diagram and attempt to 
estimate its size. In section [51 we discuss the prospects of searching non-leptonic charmless 
B c decays at LHCb and we conclude in section [HJ Two appendices are devoted to relating 
our work to results from factorisation approaches. 

2 Non-leptonic charmless B c decays as pure annihilation pro- 
cesses 

The diagram for the non-leptonic charmless B c decays is shown in fig. [T] (the case of singlet 
states will be discussed below). The initial b and c quarks annihilate into u and d or s 
quarks, which form two light mesons by hadronising with a pair of qq (q = u, d, s) emitted 
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Figure 1: Generic diagram for the non-leptonic charmless B c decays 



from a gluon. There are 32 decay channels of this kind if we consider only the lightest 
pseudoscalar and vector mesons. In the case of two outgoing pseudoscalar mesons (PP), 
there are 4 modes with strangeness one: 

K + tt°, K + r], K + T]', K°ir + 

and 4 modes with strangeness zero: 

7T+7T , 7T+7/, 7T + t/, K + K° 

The same applies for two vectors (VV) up to the obvious changes: 

K* + p°, K* + <j), K* + uj, K*°p + , p + p°, p + <j), p+u, K* + K*° 

In the case of VP decays, one can get two decay modes from one in PP decays, depending 
on the pseudoscalar meson which is turned into a vector one, yielding 8 strange decay 
modes (AS = 1 processes): 

K* + 7T°, K* + V , K* + r]', K*% + , p°K + , <f)K + , uK + , p + K° 

and 8 non-strange decay modes (AS = processes): 

P + 7T°, P + T], p+T]', K* + K°, p°7T + , 07T+, U71+ , K*°K + 

In the next section, we describe these decay channels in terms of a few reduced ampli- 
tudes using SU(3) flavour symmetry. Similar expressions have been obtained for the 
charmless B u ^ decays, which have been very useful to disentangle the rather complicated 
decay amplitudes of these decays containing many different contributions (tree, penguin, 
emission, annihilation, etc. . . ) [20, EH E2] • Comparing to the case of B u ^ decays, the 
577(3) relations for the B c decays are extremely simple, as it comes from only a single 
tree-annihilation diagram as mentioned earlier. 

The theoretical computation of the process shown in Fig. [T] amounts to determining 
the matrix element: 

(h x h 2 \H eS \B c ) (1) 
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where H. e g is the effective Hamiltonian which we discuss later on. This matrix element 
contains contributions coming from the qq state i) produced perturbatively from one- 
gaum exchange linking the dot and one of the crosses in Fig. [1]) and ii) produced through 
strong interaction in the non-perturbative regime. Which type of these two contributions 
dominates this matrix element is an important issue in the theoretical computation of the 
hadronic B decays. In many approaches to non-leptonic B decays [231 12H ESI ESI EZ] , it 
has been pointed out that annihilation diagrams may be sizable, with a large imaginary 
part, so that they have an important impact on the phenomenology of CP violation 
in B decays. Indeed, their contributions seem to be needed to bring agreement between 
theoretical computations and experimental results. There might be a significant difference 
in the annihilation contributions for B and B c decays since the B c is likely to be considered 
as a heavy-heavy system rather than a heavy-light one. We will discuss the theoretical 
estimation of the annihilation diagram in more detail in section HI 



3 Relations from SU(3) flavour symmetry 

In this section, we derive relations among the decay channels relying on the SU(3) flavour 
symmetry between u-, d- and s-quarks. Following [28J, we first write down the charmless 
B c decays in terms of the reduced amplitudes using the Wigner-Eckart theorem. 

Let us first see the possible SU(3) representation of the external states. The initial 
state, B+ is a singlet under SU(3), whereas the outgoing state is the product of two 
mesons, which can be either both in the octet representation or in one octet and one 
singlet representations. We have therefore outgoing states which transform as 

8 x 8 = 1 + 8 S + 8 A + 10 + 10* 1x8 = 8/ (2) 

where the subscripts allow one to distinguish between the three different octet representa- 
tions involved. We sandwich the operators induced by the weak interaction Hamiltonian 
between these external states to obtain the amplitudes for the B c decays. The weak 
Hamiltonian for such transitions is given by: 

n cS = [v ud v: b o AS =° + v us v: b o AS=1 ] (3) 

where the operators are: 

O AS=0 = ^(1-^)^(1-75)6 (4) 
O as=1 = ^(1-75)5^(1-75)6 (5) 

These two operators are both SU(3) octets and have the following SU(3) tensor structures: 

O AS=0 : (y,/,/ 3 ) = (0,1,1) (6) 

O as=1 : (Y, 7, 7 S ) = (0, 1/2, 1/2) (7) 
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where (Y, J, J 3 ) denotes hypercharge, isospin and isospin projection respectively. Since 
B c charmless decays involve only operators in an octet representation, one can use the 
Wigner-Eckart theorem to express all the decay amplitudes in terms of three reduced 
matrix elements: 

• a reduced amplitude S = (S^l |0 8 | |1) from the symmetric product of the two incom- 
ing octet mesons. 

• a reduced amplitude A = (8^ 1 1 C s 1 1 1) from the antisymmetric product of the same 
representations 

• a reduced amplitude / = (8/||(9 8 ||l) from the product of an octet and a singlet 
meson. 

The operator O can be O AS=0 or O as=1 . Note that the values of the reduced quantities 
S, A, I are in principle different for the PP, VP or VV final states. The Wigner-Eckart 
theorem requires one to compute the Clebsch-Gordan coefficients describing the projection 
of a given 8x8 and 8x1 final state onto the two octet operators of interest. These 
coefficients can be easily determined by combining the usual SU(2) Clebsch coefficients 
with the so-called isoscalar coefficients given in ref. [25] . 

Finally, we must consider the different symmetry properties of the out-going states 
(P and V) as discussed in ref. [30]. For PP decays, where the wave function of the final 
state is symmetric, only S contributes, apart from the case of final states containing rj 
or 7]' where both S and I are present. For VP decays, the amplitude gets contributions 
from 5" and A (and I for final states containing r],r]',u or 0). For VV decays, there 
are three amplitudes corresponding to the three possible polarisations (or equivalently 
partial waves) allowed for the outgoing state. The wave function is symmetric for S and 
D waves and antisymmetric for P wave, so that the matrix element S contributes to S 
and D waves, whereas A contributes to P waves. I contributes only to S and D waves of 
outgoing states containing <fr, uj mesons. 

A comment is in order on the mixing of the mesons containing S77(3) singlet states. 
The rj, r)', ui, <f) mesons are mixtures of the 577(3) octet (77 s or u 8 ) and singlet (77° or u°) 
flavour states 

\ V (uj)} = cos# p(l)) |77V 8 )> +sin^, ) |7/V)> (8) 
|t/(0)) = -sin^ (l ,)|r/ 8 (a; 8 )) + cos^)|r/ (^ )) (9) 

where |?7 8 ) and |u; 8 ) have the flavour composition \uu + dd — 2ss)/v / 6, and |?7 ) and \u°) 
are \uu + dd + ss) /y/3. The determination of the mixing angles 6 P>V is an important 
phenomenological issue in understanding the nature of these particles. Since we do not 
aim at a high accuracy in our SU (3) analysis, we will adopt the following values for the 
mixing angles which are not very far from the phenomenological determinations: 

tan# p = — —, tan 6^ = \pl. (10) 
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These angles correspond to the ideal mixing for the vector sector: 

uj = (uu + dd) I V2 <f) = ss (11) 

and also yield a simple expression of the pseudoscalar mesons: 

rj = (uu + dd- ss)/V3 rf = (uu + dd + 2ss) / 'y/E (12) 

The non-ideal mixing of the (r],r)') mesons is linked to the U(1)a anomaly (see, e.g., 
refs [3"T| E2] ) • This value of the pseudoscalar mixing ang le e p = arctan(2 v / 2)" 1 ~ -19.5° 
is close to phenomenological determination, e.g. from the J/tp radiative decays, 9 P ~ 
—22° [33]. Let us stress that for the light mesons, we take the same phase conventions 
as in ref. [22], so that some amplitudes have a minus sign with respect to those obtained 
from ref. [22] 0. 



3.1 PP modes 

Taking into account the Clebsch-Gordan coefficients together with the issue of octet- 
singlet mixing, we obtain the following amplitudes for the PP modes, we have 

Amplitude 


3 qPP 

4 gPP _|_ V2jPP 
3\/5 3 
-l^jlS p P _|_ ijPP 

Here and in the following tables, these amplitudes must be multiplied by Gp/y/2 and also 
the appropriate CKM factor V u r>V* b with D = d or s. We notice the relations 

A(B+ -> A'°7r+) = V2A(B+ -> A+7T ) = XA(B+ -> K+K°) (13) 

with the Cabibbo-suppressing factor A = V us /V u d- The above relations are valid in the 
exact SU(3) limit (for instance, we have S pp = S K+7T ° = S K ° n+ = S K+K ). Obviously, 
these relations have some interest only if the size of SU(3) breaking remains limited - we 
will discuss this issue in Sec. HI 



Mode 


Amplitude 


Mode 


A+7T 


Hi" qPP 

V 10 


71+71° 


K°TT + 




K + K° 


K+T] 


2 qPP 1 V2 tPP 


7i + r] 


K+T]' 


1 qpp _l 4 tpp 


7i +r q' 



*In detail, (—u, d, s) transform as an anti-triplet [3U], which means that there is a (—1) phase between 
the conventions of refs. [55] and [H] for the pseudoscalar mesons ir~ , ir°, K ~, r}, rf and the vector mesons 
p~ , p°, K*~ , (j>, u>. We have multiplied all the amplitudes by a further (—1) factor, so that the differences 
between our results and those obtained using ref. [29] are limited to a (—1) factor for the decay amplitudes 
for K°ir + , K + K° and their vector counterparts. 
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3.2 VP modes 

For the VP modes, we have for the strange modes 



Mode 


Amplitude 


Mode 




1 fscVP, 1 A VP 


uK + 


p°K + 


1 fzqVP 1 aVP 

2V5 5 27/3^ 


(f)K + 


K*+7] 


— ly/I^VP _|_ V2j^VP _j_ IjVP 


p + K° 


K *+ v > 


_J_gVP _ lj{VP _|_ 2V2 jVP 
6^/5 6 3 





and for the non-strange modes 



Mode 


Amplitude 


Mode 


p+7T° 


1 AVP 


C<J7T + 


p°7T+ 


~M VP 


4>7T + 


p+?7 




K*+K° 


p+7?' 


l_gVP _|_ gygjVP 

3\/5 3 


K*°K+ 



providing the simple relations 



A(B- 



K*\ + ) 
♦ P + K°) 



Amplitude 



1 



5 



vp !—A VP + -I VP 

2^/3 y 3 



3 r-yp 1 aVP 
3 QVP i 1 4 VP 

To* 5 + 7e 



Amplitude 



3 cFP 
10 



3 qVP 1 1 a VP 



iT+TT ) = AA(^ 

p°K+) = \A{Bt 



K* + K°) 



(14) 
(15) 



It should be noted that the amplitude J yp can be significantly different for the pro- 
cesses involving the vector singlet (0, u>) and the pseudoscalar singlet (77, 7/) since it is 
known that the latter should receive a contribution from the anomaly diagram. This 
could induce a significant breaking of the above relations for channels involving rj, rj'. 



3.3 VV modes 

For the VV modes, we have three different configurations for the outgoing mesons, labeled 
by their (common) helicity. The left-handedness of weak interactions and the fact that 
QCD conserves helicity at high energies suggest that the longitudinal amplitude should 
dominate over the transverse ones (corresponding to helicities equal to ±1). These helicity 
amplitudes can be combined linearly into S, P and D wave amplitudes, and in particular, 
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the longitudinal amplitude is a linear combination of only S and D waves. 



Mode S, D Amplitudes P Amplitude 

K*+n° flLqvv J-A VV 

A P V 10^0,2 

&-*+, , i qvv I 2 rvy i aw 

n~qvv_i_ Rjvv n a vv 

A ^ yT5 D 0,2 + y 3 J 0,2 y 3 71 1 

K-*o„+ Ilsqvv i 4W 

A P V 5^0,2 75^1 

p+p° y^A^ 

„+, , FVqvv I 2 ryy n 

P w y T5'- , 0,2 + 75^0,2 U 
2_qVV , RjVV n 



/15^0,2 ~ Y 3 0,2 

zV*+r7*0 IHqVV 1 

A A V 5^0,2 

where the subscript denote the partial wave under consideration I = 0, 1, 2. In particular, 
we have the interesting relations 

A(B+->K*°p + ) = V2A(B+ -»• K* + p°) (16) 
AA(5+ K*+K°) = V2(-l) e A(B+ ^ K* + p°) (17) 



3.4 Zweig rule 

In the above expressions, the normalization between S, A and / amplitudes are different: 
the former is related to the 8x8 representation and the latter to 1 x 8. One may relate 
these two amplitudes by means of the Zweig rule for the AS* = processes involving 0. 
At the level of quark diagrams, one can see that the B^t — > 07r + (p + ) process cannot come 
from the diagram in fig. [T]- since is a pure ss state. The only production process come 
from non-planar diagrams diagram similar to fig. [U, but with a different combination of 
quarks into the outgoing mesons: the ud quarks produced from the W go into 7r + (p + ) 
whereas the meson is made of a ss-pair produced from vacuum. Such a non-planar 
diagram is expected to be suppressed, especially for vector mesons such as (at least 
three gluons are needed perturbatively, and it is 1/N C suppressed in the limit of a large 
number of colours). 

Assuming that the amplitudes for B+ — > 07r + and B+ — > 0p + vanish, one obtain the 
following relations: 

I VP = \fls VP IU = ]fl^ (18) 



S 



providing simpler expressions for the following VP decay amplitudes 



Mode 
K*+r) 

K*+7]' 

luK+ 
<f)K + 



Amplitude 



Mode Amplitude 



a/I ^4 VP 



3 qVP 

175 



3 gVP 



3 cVP 
10° 



1 A VP 

J_AVP 

V6 A 



P + V 
p + rj' 



LU7T 



(f)TT~ 



'2gVP 
l gyp 

3gVP 



so that 



A(B+ -> p+ri) = V2A(B+ -> p+r/) (19) 

Although we assume here that the I VP amplitude is the same for the processes involving 
the vector singlet (<f),w) and the pseudoscalar singlet (77,77'), as required by 577(3) sym- 
metry, this assumption is broken for the pseudoscalar singlets due to the anomaly (seen 
for instance in the mass difference between r/ and 77'). At the quark level, the detached 
diagram for the pseudoscalar singlet states (77, 77') cannot be neglected since there is the 
well-known anomaly contribution modifying the previous relation: 



V2 



A{Bj 



P V 



2V2 



AI 



VP 



(20) 



where AI VP denotes a potentially large l/iV c -suppressed anomaly contribution. 
Keeping the same caveat in mind, we can simplify some VV amplitudes 

Mode S, D Amplitudes P Amplitude 





/~3qVV 
V 10°0,2 


1 AVV 


K*+(j) 


[3 qVV 
V 5^0,2 


[2 aVV 
V 3 1 


p + u 


f&qVV 
V 5 D 0,2 





p+0 









with the obvious relations among partial waves. 



4 Estimating the branching ratios 

As mentioned in section [21 a precise estimate of the matrix element for the annihilation 
diagram is an important theoretical issue in B physics. Although the domination of the 
one-gluon exchange diagram has been argued in various theoretical frameworks for B u ^ s 
decays [251 G3], it has not been investigated whether one-gluon exchange or other (non- 
perturbative) contributions dominate in B c decays. In this section, we provide branching 



9 



ratio estimates for the non-leptonic charmless B c decays in two ways, by using experimen- 
tal data on pure annihilation B decays and by relying on the one-gluon picture a la QCD 
factorisation. The branching ratios can then be readily obtained by the usual formulae 

J [M 2 B - (mi + m 2 ) 2 pi| - (m x - m 2 ) 2 ] 
Br(B c - hiht) = * rtotl67rM 3 C \{hM^\B c )\ 2 (21) 

and the expression obtained in the section El in terms of the reduced amplitudes is related 
to the matrix element through 

{(h^lTUsWl = G F /V2\V ud{s) V; b \ x \R(h h h 2 )\ (22) 

where the reduced amplitude R(hi, h 2 ) involves the amplitudes listed in sec. El expressed 
in terms of (S, A, I). 



4.1 Estimate from annihilation process 

There are two pure annihilation processes observed in heavy-light B decays [3U ESI EH] : 

Br(B° — > K + K~) = (0.15+°;lo) x lO" 6 (23) 
Br(B° -> DJK+) = (3.9 ± 2.2) x 1(T 5 (24) 

Although the large experimental errors do not allow us to draw any firm conclusion, 
these data seem to indicate that the annihilation contribution is not negligible. We may 
attempt to use these decay channels to very roughly estimate the size of the non-leptonic 
charmless B c decays. Since we are interested in charmless final states, let us compare 
the B c -> K + K and the B° -> K+K~ processes. Assuming naive factorization between 
initial and final states, the final-state contribution cancels out when taking the ratio of 
the amplitudes. As a result, we find: 

Br(B c -> K+K°) 

The factor ^ represents the difference due to the fact that the B° K + K~ process comes 
from a diagram similar to Fig. [1] but the W boson propagates in the t-channefj]. In the 
one-gluon picture, £ = C\jC<i ~ 4 whereas it might be smaller once non-perturbative 
effects are included. Indeed these effects would add new contributions to both decays, 
but the B c decay would be more affected relatively, since its Wilson coefficient in the one- 
gluon picture is smaller. This very naive argument leads to the relation between these 
two branching ratios: 

1 o ss 1 

Br(B c -> K + K°) ~ Br(B° -> K + K~) x — — > Br(B° -> K + K~) x 7.5 (26) 




^Here we neglect the small penguin diagram contribution to the B° — > K + K decay. 
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The estimate for Br{B c — > K + K ) using this relation depends on the result on Br(B° — > 
K + K~), which should be improved in the near future. Taking the current central value 
of Br(B° -> K+K~), we find a lower limit of Br(B c -> K+K°) at the order of KT 6 . 
Using this result, one can estimate the Wigner-Eckart reduced matrix elements S pp : 

S pp > 0.085 GeV 3 [(PP) = (K°K + )} (27) 

where we used the following CKM central values: V u b = 0.0035, V c b = 0.041 [37] . 

As mentioned before, there is no good reason to assume that S pp , S PV and S vv 
should be related. Since we are only looking for order of magnitudes, we will assume as a 
dimensional estimate that \S PP \ ~ \/2\S PV \ ~ (S^l- We emphasise that these relations 
have no strong theoretical supports and are just meant as a way to extract order of 
magnitudes for the branching ratios. In addition, we assume the Zweig rule to determine 
the singlet contributions I and we neglect the antisymmetric contributions A, as well as 
transverse VV amplitudes. This provides the following branching ratios of interest, for 
instance: 

[B d annihil] BR(B C -> (j)K + ) ~ C(10~ 7 - 10~ 8 ), BR(B C -> K*°K + ) ~ C>(10~ 6 ) 

BR{B C -> K°K + ) ~ C(10~ 6 ), BR(B C -> K*°K* + ) ~ O(10" 6 ) (28) 

The suppression of the <j)K + channel is due to the small CKM factor for the AS = 1 
processes, Cabibbo-suppressed compared to AS = 0. 



4.2 Estimate from one-gluon exchange model 

A second method consists in a model based on one-gluon exchange, in close relation with 
the model proposed in QCD factorisation to estimate annihilation contributions for the 
decays of heavy-light mesons. In this method, described in more detail in App. A, the 
matrix element in eq. ([T]) can be given as: 

(h 1 h 2 \'Hrt\B c ) = t^V; b V ud{s) N hlh2 b 2 (hu h 2 ) (29) 

where 

N hlh2 = f B J h J h2 hih, h 2 ) = ^LCiAKhh,). (30) 

The function A\(h\h2) is estimated as the convolution of the kernel given by one-gluon 
exchange diagrams and the distribution amplitudes of the initial and final state mesons. 
While the detailed computation of this function can be found in Appendix A, we would 
like to emphasize a few differences in B c decays comparing to the B u d decays that we 
found; i) the B c decays are much simpler, since the only operator contributing is 2 , 
and there is only one combination of CKM factors (V* b V U D where D = d,s depending on 
the strangeness of the outgoing state), ii) the long-distance divergences, which prevents 
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us from estimating the annihilation contribution in B u d decays, do not appear in B c 
annihilation. 

We can give the numerical results for the following channels: B c — > <f)K + ,B c — > 
K*°K+,B C -> K°K + ,B C -> K°*K* + . We start with the function b 2 (h 1 h 2 ) which turns 
out to be quite SU(3) invariant: 

6 2 (0^+) = 1.5, b 2 {K*°K + ) = 1.4, (31) 
b 2 (K°K + ) = 1.6, b 2 (K°*K* + ) = 1.0 (32) 

Notice that 6 2 (VP) = &2CPV), and thus A yp = 0. One can see that the S77(3) breaking 
is rather small as argued in Appendix A, while the difference between PP and VP(VV) 
modes can be as large as +13(+38)%. We next list the normalization factors: 

N+ K+ = 0.014 GeV 3 , N K .o K + = 0.014 GeV 3 , (33) 
N^o K+ = 0.010 GeV 3 , N K°* K *+ = °- 019 GeV3 ( 34 ) 

For these particular decay channels (no exchange between ir and K), the £77(3) breaking 
is also small while the difference between PP and VPiVV) modes can be as large as 
-40(-90)%. Our numerical value for £, I, A amplitudes for the above processes are: 

S VP = 0.036 GeV 3 , [VP = (f)K + ,K* K + } (35) 
S pp = 0.021 GeV 3 , [PP = K°K + ] (36) 
S% v = 0.025 GeV 3 [VV = K°*K* + ] (37) 

where we have neglected transverse VV amplitudes. We present the relation between 
Nhxhi and b 2 (hih 2 ) functions and the Wigner-Eckart reduced matrix elements S, I, A in 
App. B. Assuming the SU(3) breaking effect is negligible in b 2 (hi, h 2 ), this relation allows 
us to estimate all the other S, I, A amplitudes with the values given in Eq. (I3"2l and the 
known values of decay constants for each final state. 

We obtain finally the following values of branching ratios: 

[One-gluon] BR(B C -> <f)K + ) = 5 x 10~ 9 , BR(B C -> K*°K + ) = 9.0 x 10~ 8 (38) 
BR(B C -> K°K + ) = 6.3 x 10^ 8 , BR(B C -> K*°K* + ) = 9.1 x 10^ 8 (39) 

The contributions to p ir + , 07r + and p + <f) vanish in our approximations, which means that 
these power-suppressed decays must have significantly smaller branching ratios than the 
above ones. We do not quote any error bars on these results on purpose: we can easily 
estimate the uncertainties coming from our hadronic inputs, but certainly not the sys- 
tematics coming from the hypothesis underlying our estimate (one-gluon approximation, 
asymptotic distribution amplitudes, neglect of l/m b and l/m c suppressed corrections, 
neglect of soft residual momentum of the heavy quarks in the B c meson). 
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4.3 Comparison of the two methods 

Our estimates of the branching ratio for e.g. the B c — > K + l{ in the above two ways are 
not consistent. This is clearly because two methods are conceptually different: 

• The method based on annihilation treats the charm quark as massless. It takes 
into account some of the non-perturbative long-distance effects expected to occur 
in Bd and B c decays, but treats in a very naive way the relation between matrix 
elements of the operators 0\ and O2. It relies also on extremely naive assumptions 
concerning the respective size of matrix elements for PP, VP and VV modes. 

• The method based a perturbative one-gluon exchange treats the charm quark as 
heavy. It assumes the dominance from a specific set of diagrams computed in a 
perturbative way, but it provides a consistent framework to perform the estimation. 

It is well known that both kinds of estimates yield rather different results. This is illus- 
trated by the fact that the estimate of B^ — > K + K~ in the annihilation models of QCD 
factorisation [25] (around 10 -8 with substantial uncertainties) and perturbative QCD [HE] 
is one order of magnitude below the current experimental average. Therefore, it is not 
surprising that our two methods yield branching ratios differing by a similar amount. 
There are well-known cases where final-state interaction can increase significantly esti- 
mates based on factorisation, for instance B — > K\ c [3"§l HO] or £)+ — > p°7r + |41j . An 
observation of the non-leptonic charmless B c decays will certainly have an important key 
to clarify such a controversy as well as further theoretical issues in computational methods 
for the annihilation diagram. 

5 Search prospect at LHCb 

The LHC pp collider with the center of mass energy of 14 TeV has a large cross section 
for the bb hadro-production, which can be followed by the production of not only B°, 
B + and B® mesons but also other b hadrons such as, Af, and B c . The subtraction of the 
production of other known b hadrons [46J 

bb -> (B d : B U :B S : A' b s) ~ (42.2 ± 0.9% : 42.2 ± 0.9% : 10.5 ± 0.9% : 9.1 ± 1.5%) (40) 

leads to bb — > (B c ) to be less than one %. The LHCb experiment [17], which is dedicated 
to B physics analyses with its optimized trigger scheme, allows one to detect the 6-decay 
modes into hadronic final states. 

The theoretical estimate of the B c cross section is still under scrutiny. For the dom- 
inant gg-fusion process, there are two possible mechanisms: gg — > bb followed by the 
fragmentation, or b — > B c bc. It is found that the latter dominates in the low-pr region 
which corresponds to the LHCb coverage [HJII9]. The 0(a^) computation of the b — > B c bc 
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process predicts a(pp — > B^X) ~ 0.3 — 0.8 yh [HU] where the error comes from the un- 
certainty on theoretical inputs such as the choice of the a s scale and the B c distribution 
function. Additional systematics could come from higher-twist and radiative corrections. 
In the following, we follow the LHCb value for the cross section cr(B c ) = 0.4 pb but it 
must be noted that this value may be affected by a large uncertainty. 

We can now estimate the expected sensitivity for a specific channel. First, let us 
discuss which channel has the best potential for the detection. The best trigger and 
reconstruction efficiencies with a large signal over background ratio can be achieved by 
the charged K- and/or 7r-tags (and by avoiding low-pr neutral particles) at LHCb. Since 
the initial B c carries an electric charge, all two-body PP final states contain one neutral 
particle. The same remark applies for the VV channels when one considers the subsequent 

decays of the vector particles into pairs of pseudoscalars. In this respect, PV channels 



such as B+ -> <pK + , K K + , K 7r + , p K + , p°7r + , (f)ir + are the best candidates using 
the vector meson decays, <p K + K~, K — > K~tt + , p° — > 7r + 7r~, leading to three 
charged tracks. Among these subsequent decays, the small widths of and K make the 
reconstructing particularly easy comparing to e.g. p°. On the theoretical side, our Zweig 
rule argument forbids the B+ — > 07r + , whereas the B+ — > p°7r + channel comes only from 
the A (asymmetric) amplitude which is also subdominant. Finally, taking into account 
the fact that the AS = 1 channels are Cabibbo suppressed, we draw the conclusion that 
the Bj — > K K + channel might be the best candidate for the detection. 

Since the selection criteria and trigger efficiencies are different for each channel, de- 
tailed simulations are necessary in order to estimate the expected sensitivity for different 
channels. For example, such a study has been done for B c — > J/tp7r + [501 EI]- From 
the expected branching ratio Br(B c — > J/tp7r + ) ~ 1 %, it was deduced that over a 
thousand of events are expected after the one year run of LHCb. By scaling this ob- 
servation to the processes of interest, we can very roughly estimate that an assumption 
of Br(B+ -> K K+) = 10~ 6 yields a few events per year at LHCb. The analysis of 
LHCb data will thus allow to set first experimental limits on the non-leptonic charmless 
B c decays, and give hints on annihilation mechanisms in these decays. 

6 Conclusions 

In this paper, we have discussed non-leptonic charmless B c decays into two light pseu- 
doscalar or vector mesons. It turns out that a single tree annihilation diagram is respon- 
sible for all 32 processes, providing an interesting testing ground for annihilation. After 
discussing general aspects of the charmless B c decays, we have shown that the very simple 
nature of these decays allows us to describe them in terms of a few reduced amplitudes 
by exploiting SU(3) flavour symmetry to relate various PP, PV and VV modes. 

In order to discuss a possible search for charmless non-leptonic B c decays at LHCb, 
we have proposed two different theoretical estimates of these reduced matrix elements, 
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either by comparison with annihilation processes or by a perturbative model based on 
the exchange of one-gluon. The two models yield a rather wide range of branching ratio 
predictions, from 10~ 6 to 1CT 7 . The LHCb experiment has the potential to observe some 
of the decay channels (such as B c — > (pK + , K K + ) if the branching ratio is at the larger 
side of these estimates. 

From the theoretical point of view, a better understanding of annihilation diagrams 
is particularly important. They are often assumed to play a significant role in decays 
of heavy-light mesons, but they occur jointly with other kinds of diagrams, making it 
difficult to assess precisely their size. Furthermore, for the theoretical estimates of B u ^ 
annihilation diagram in the QCD factorisation, there is an additional uncertainty caused 
by the infrared divergence occurring in its computation. It is worth mentioning that we 
found that such a divergence does not occur in the case of the B c annihilation diagram 
suggesting that predictions from models a la QCD factorisation for the B c decays should 
be more precise, and thus easier to confirm or reject. 

On the other hand, it has been discussed that the annihilation diagrams may be 
enhanced by long-distance effects such as final-state interactions. Although only limited 
models of such effects have been proposed either for D or for B decays (the former likely 
more affected than the latter by such enhancements) [HI 1391 14*0] . the observation of an 
unexpectedly large branching ratio for the B c annihilation would call for a reassessment 
of such long-distance contributions. An observation of charmless non-leptonic B c decays 
at LHCb will certainly provide substantial information on these models, in complement 
with the observation of other decays such as Ba — > K + K~ or B s — > 7r + 7r~. 
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A Short-distance model for weak annihilation 

As highlighted in the introduction, weak annihilation plays a significant role in B u ^,s 
non-leptonic decays, but it is difficult to estimate it accurately. A model to estimate this 
contribution was provided in the framework of QCD factorisation [4*21 125]. relying on the 
following hypothesis : 

• the diagrams are dominated by the exchange of a single gluon, whose off-shellness 
is typically of order 0(^fkrr^ } ) 
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• hadronisation effects are taken into account through light-cone distribution ampli- 
tudes (generally taken in their asymptotic form) 

• soft components are neglected. 

Being power-suppressed in the heavy-quark limit, the weak-annihilation contributions 
to B U) d )S non-leptonic decays cannot be factorised in short- and long-distance effects (in 
general). Their evaluation within this rough model exhibits endpoint divergences, which 
signals the presence of long- distance contributions not taken into account properly. The 
divergent integrals were regularised on the basis of dimensional analysis, which induces a 
significant uncertainty on the estimate of the annihilation contribution. 

We can follow a similar method to estimate annihilation in the case of the B c decay. 
Concerning the B c meson, we work in the limit where both b and c quarks are heavy 
(keeping m c / nib fixed) and we set the momentum of the valence quarks to p^ = m b v ^ and 
p^ = rricV 11 , neglecting the soft components of the heavy-quark momenta (and consistently 
setting Mb c = m c +nib). Since we neglect the soft components of pb andp c , the integration 
over the B c meson distribution amplitude is trivial and yields fs c - The diagrams to 
compute are not very difficult and correspond to a gluon emitted from the b anti-quark or 
the c quark from the B c meson and converted into a light quark-anti-quark pair. Following 
ref. [25], we find in the case where (Mi, M 2 ) = (P, P), (P, V), {V, V) : 



A\(M 1 M 2 ) = 7ta s I dxdy (41) 

1 1 



<t>Ah(y)<t>M 3 (x) 



+r Ml r M ^ mi (y)4> m2 (x) 



_y[(x + y)z b -xy] x[(x + y)z c - xy] 
2(1 - z b ) 2(1 - z c ) 



(x + y)z b -xy (x + y)z c - xy 



If (Mi, M 2 ) = (V, P), one has to change the sign of the second (twist-4) term above. 4>m 
and 4> m are twist-2 and twist-3 two-particle distribution amplitudes of the meson M, and 
r M is the normalisation of the twist-3 distribution amplitude. In the case of pseudoscalar 
mesons, we have: 

r* = 2ml r* = , 2ml , (42) 

nib x 2m q nib(m g + m s ) 

responsible for the chiral enhancement of twist-4 contributions for pion and kaon outgoing 
states. In the case of vector mesons, we have: 

r v = 2 ^ f ± (43) 
ni b Jv 

z b and z c denote the relative size of the b and c-quark masses: 

z h = ; z c = 1 - z b = ■ (44) 

m b + m c m b + m c 



1(3 



Their appearance allows one to distinguish the diagram of origin (corresponding to a 
gluon emitted from the b or the c quark line). 

Eq. (1411) is in agreement with the expressions obtained in refs. [25] and [13] in the 
limit Zb — > 1 (and z c — > 0). To simplify further the discussion, we take the asymptotic 
expression for the distribution amplitudes 

(f) P (x) = 6x(l-x) <f) V {x) = 6x(l -x) (j> p (x) = l <p v {x) = 3{2x - 1) (45) 

The structure of the singularities in the kernel is due to the propagator of the gluon and 
seems quite complicated. But if we take as an example 



i 1 j 



dx dy— - — — dx r log 

(x + y)z — xy J z — x 



z — X + xz 



xz 



(46) 



The function to be integrated is continuous at x = z and has integrable singularities for 
x = and x = z/(l — z). The integration can therefore be performed without problem as 
long as z is different from 0, 1/2 and 1 (no coalescence of singularities). For the B c meson, 
it means that the twist-2 and the twist-4 contributions have no endpoint singularities and 
yield finite integrals. Therefore, there is no need to introduce models to regularise the 
divergent integrals like in the case of heavy-light mesons. 

The corresponding expressions for the four cases are slightly tedious, but they can 
be approximated to a very good accuracy through low-order polynomials in 5, where 
Zb = 0.76 + 5 and z c = 0.24 — 5 (corresponding to mj = 4.2 and m c = 1.3 for 5 = 0): 

A\(PP) = 7ra s [(-22.83 + 4.845 + 808.35 2 + 25075 3 + 34255 4 ) (47) 

+r Ml r M2 (-8.23 - 3.655 + 73.65 2 - 16. 15 3 + 35755 4 + 160075 5 )] 
A\(PV) = 7ra s [(-22.83 + 4.845 + 808.35 2 + 25075 3 + 34255 4 ) (48) 

+r Ml r M2 (-19.15 - 123.55 - 1305 2 + 58.75 3 + 79825 4 + 397785 s )] 
A\(VV) = 7ra s [(-22.83 + 4.845 + 808.35 2 + 25075 3 + 34255 4 ) (49) 
-r Ml r M2 (2.44 + 222.45 + 15655 2 + 33865 s + 58245 4 - 808315 5 - 4219275 6 )] 



Within the set of approximations performed here, A\(PV) and A\(PV) are identical. 

We use the above formulae to estimate a few branching ratios. We take our inputs for 
the vector decay constants and the Wilson coefficient C 2 (\/miAh) = —0.288 from ref. [4*4] . 
and the rest of our inputs from ref. [23] . We take the value of the B c meson decay constant 
fBc = 395 MeV taking the central value from ref. [3] . 

Let us comment on the SU (3) breaking, which can be included in this QCD computa- 
tion. We do not have the 5 , '7(3)-breaking effects coming from the distribution amplitude: 
for instance, a small m s correction makes the K and distribution amplitudes slightly 
asymmetric. On the other hand, we have the breaking effect in the chiral enhancement 
parameter r M . The SU(3) breaking (e.g. comparison of r n ,r K or r p ,r K * ,r^) turns out 
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to be relatively small. Another SU(3) breaking arises from the following decay constants 
in the normalization factor N^hi [?> 

f n = (130.4 ± 0.2)MeV, f K = (155.5 ± 0.8)MeV (50) 
f p = (216 ± 3)MeV, f K . = (220 ± 5)MeV, f u = (187 ± 5)MeV, fa = (215 ± 5)MeV. 

The SU(3) breaking in the decay constants for the vector mesons is rather small while 
there is a 24% difference in ir and K decay constants. 



B Identification with results from QCD factorisation 

The expressions for all the decay channels considered in Sec. [3] can be recovered from 
refs. [25] , [13] and [15] if we identify between the Wigner-Eckart reduced matrix elements 
S, I, A and the 2 reduced coefficients b 2 . In these references, one must take the expres- 
sions for the decay amplitudes of B u decays into the relevant final state, and pick up the 
2 contribution, which is the only remaining one once the related B c decay is considered. 
If we perform this identification, we obtain: 

S PP = ^N PP b 2 (PP) (51) 

I PP = ^N PP b 2 (PP) (52) 

S VP = ^N VP (b 2 (PV) + b 2 (VP)) (53) 

A VP = ^N VP (b 2 (PV)-b 2 (VP)) (54) 

I VP = ^N VP (b 2 (PV) + b 2 (VP)) (55) 

iVV _ . / 5 at uS,D, 



S V S,D = ^-N VV b^(VV) (56) 

A V P V = (57) 
I V s V d = \J\Nvvb S 2 D {VV) (58) 
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